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We report the latest results of the ZEUS collaboration on the photoproduction of D* mesons in a low W 
range. Differential cross sections as function of and rj D are measured and compared with several recent NLO 
pQCD calculations. The differential cross-sections in a restricted kinematical region are higher than the NLO 
calculations, in particular in the forward (proton) direction. A recent pQCD model (BKL) describes the data 
reasonably well. 



1. INTRODUCTION 

Heavy quark photoproduction can be used to 
probe pQCD calculations with a hard scale given 
by the heavy quark mass and the high trans- 
verse momentum of the produced parton (toq 3> 
Aqcd)- Two types of NLO calculations with 
different approaches are available for comparison 
with measurements of charm photoproduction at 
HERA. The massive charm approach [ [jj assumes 
light quarks to be the only active flavours within 
the structure functions of the proton and the pho- 
ton, while the massless charm approach [||, |§ also 
treats charm as an active flavour and is thus only 
valid for p± 3> m c . 

The data taken by the ZEUS collaboration dur- 
ing 1996/1997 corresponds to an integrated lumi- 
nosity of about 37 pb -1 . In a subsample of about 
17 pb -1 a small calorimeter positioned along the 
beam pipe was used to tag low W events, 80 < 
W-y P < 120 GeV. The results of the high W region 
(130 < W JP < 280 GeV) have been published be- 
fore [ U and will not be shown here. This is the 
first presentation of our low W results. Charm 
was identified by the observation of I?* (2010) 
mesons, which were reconstructed in the follow- 
ing decay modes: D* + — ► D°irf — > (K~7r + )nf 
(Br = 0.0262 ± 0.0010) and D*+ -> D°ir+ -> 
(A^7r+7r+7r-)7r+ (Br = 0.051±0.003) and charge 
conjugates. The kinematic range studied was 



pY" > 2 GeV and -1.5 < if* < 1.5 for the 
high W region, and 2 < p^ < 8 GeV and 
-1.0 < 1]°" < 1.5 for the low W region. The 
pseudorapidity is r] D = — ln(tan|), where 8 is 
the polar angle with respect to the proton beam 
direction. 

Charged tracks were measured in the central 
tracking detector. Cross sections were calculated 
in the photoproduction range of photon virtual- 
ities Q 2 < lGeV 2 ( Q 2 < 0.015 GeV 2 for the 
tagged data). 

2. D* RECONSTRUCTION 

D* events have been selected by means of the 
mass difference (AM) method [ |5j. In the high 
W region we have observed [ £§ 3702 ±136 D* 's in 
the D° -> (Avr) decay mode with pf * > 2 GeV, 
and 1397 ± 108 in the (Ktttttt) decay mode with 
pJf > 4 GeV (M(D°) = 1.80-1.92 GeV). In the 
low W region we triggered only the (Kir) decay 
mode, and observed 550 ± 36 D* events in the 
range 2 < ' < 8 GeV (Fig. 1). All tracks were 
assumed to be pions and kaons in turn; wrong 
charge D* combinations [ [I| were used as a back- 
ground distribution (dashed curve in Fig. 1), nor- 
malized outside the signal region. 

3. D* CROSS SECTIONS AND COM- 
PARISON WITH CALCULATIONS 
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The D* differential cross sections in the low W 
region are shown in Fig. 2 to 5. For comparing the 
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Figure 1. AM distribution for the events inside 
the D° mass region. The wrong charge distribu- 
tion is shown as a dashed histogram. 
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experimental data to the NLO QCD calculations, 
we have used the D* branching value measured 
by OPAL [§ f(c -> D*+ + ...) = 0.222 ± 0.014 ± 
0.014. For the charm fragmentation to D* the 
Peterson [ |7j fragmentation function was used: 

n < \ at z(l~z) 2 p D , 
Dc{z) = N [(l-zr + ez]i> Z = — - 

In the massive calculation e = 0.036 was ob- 
tained from a recent fit of Nason and Oleari [ 
H to ARGUS data. Alternatively, the Peterson 
fragmentation was replaced by fragmentation ef- 
fects estimated by a leading order Monte Carlo 
(Pythia). Initial and final state radiation were 
not included. The results of both calculations for 
the low W region are shown in Figs. 2 and 3. 
The cross sections are compared with NLO QCD 
massive calculations using MRSG for the proton 
structure function (SF) and GRV-G HO for the 
photon. The theoretical massive calculations are 
below the data, in particular in the forward (pro- 
ton) direction, although the Pythia fragmenta- 
tion slightly improves the agreement. 

A comparison with massless calculations [ |2| , 
which are expected to become valid mainly at 
higher p® , is shown in Fig. 4 for several pho- 
ton structure functions. Some sensitivity to the 
photon SF seems to be present, but the excess in 



Figure 2. ZEUS tagged data compared to massive 
NLO calculations of da /dp® . 



the forward direction is evident. The structure 
function GS-G-96 HO describes the data best. 

Recently Berezhnoy, Kiselev and Likhodcd 
(BKL) have suggested a new model for describ- 
ing D* photoproduction [ ||. In this tree level 
pQCD O(ao i s ) calculation, they hadronize the 
(c, q) state produced in pQCD, taking into ac- 
count higher twist terms at p± ~ m c . Thus the 
model is supposed to be valid over the whole p± 
range studied. No explicit resolved component 
is used. Singlet and octet color states both con- 
tribute to D* production. The color state ra- 
tio 0(8)/0(l) is a free parameter in this model 
and was tuned to the ZEUS untagged results [ ||] , 
yielding a value of 1.3. 

Comparison of these calculations, for the same 
Octet/Singlet mixture, with the ZEUS tagged low 
W data is shown in Fig. 5. A better agreement 
with the data is observed than that for the NLO 
calculations. 
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: :: : :: PYTHIA fragm. 
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Figure 3. ZEUS differential cross sections da/dr/ 
compared to the massive NLO predictions. The 
dashed/dashed-dotted lines correspond to the 
Pythia fragmentation. Parameter sets are as in 
Fig.2. 
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Figure 4. ZEUS data and massless NLO predic- 
tions [ U of differential cross sections for several 
photon structure functions. Peterson fragmenta- 
tion was used with e = 0.116. 
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Figure 5. Comparison of the ZEUS low W data 
with the BKL model. The Octet/Singlet ratio is 
1.3, as tuned for the high W region. 



